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Abstract 

We consider the heat equation in a straight strip, subject to a combi- 
nation of Dirichlet and Neumann boundary conditions. We show that a 
switch of the respective boundary conditions leads to an improvement of 
the decay rate of the heat semigroup of the order of t~ x ^ 2 . The proof em- 
ploys similarity variables that lead to a non- autonomous parabolic equa- 
tion in a thin strip contracting to the real line, that can be analyzed on 
weighted Sobolev spaces in which the operators under consideration have 
discrete spectra. A careful analysis of its asymptotic behaviour shows 
that an added Dirichlet boundary condition emerges asymptotically at 
the switching point, breaking the real line in two half- lines, which leads 
asymptotically to the 1/2 gain on the spectral lower bound, and the f~ 1//2 
gain on the decay rate in the original physical variables. 

This result is an adaptation to the case of strips with twisted boundary 
conditions of previous results by the authors on geometrically twisted 
Dirichlet tubes. 
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1 Introduction 

We consider the heat equation 

Ut - Au = (1) 

in an infinite planar strip f2 := R x (—a, a) of half-width a > 0, subject to 

Dirichlet boundary conditions on V® := (-co, 0) x {—a} U (0, +oo) x {a} , 
Neumann boundary conditions on := (0, +oo) x {—a} U (— oo, 0) x {a} , 

and to the initial condition 

4,0) = « ei 2 (!]). (2) 

This model is considered as a 'twisted' counterpart of the explicitly solvable 
problem given by (see Figure [1]) : 

Dirichlet boundary conditions on Tq := (-co, +oo) x {—a} , 
Neumann boundary conditions on Tq := (— oo, +oo) x {a} . 

Henceforth we shall use the common subscript 

e e {o,tt} 

when we want to deal with the two problems simultaneously (the value of 9 
suggests the rotation angle giving rise to twisting/untwisting). 



Figure 1: Planar strips with untwisted (left) and twisted (right) boundary con- 
ditions; the thick and thin lines correspond to Dirichlet and Neumann boundary 
conditions, respectively. 

The solution to ([I])-© is given by u(t) = e At>t u , where e Aet is the semi- 
group operator on L 2 (f2) associated with the Laplacian — Ag determined by the 
respective boundary conditions (depending on 9). 

The operators — and — Ao have the same spectrum 

o-(-A e ) = er css (-A ) = [Ei,oo) , where E x := (-^V . (3) 



Consequently, for alH > 

lL 2 (0)->L 2 (n) 



V4a 



p A 9*|| — p-Elt (A\ 



irrespectively of the value of 9. 
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In this paper, we are interested in additional time decay properties of the 
heat semigroup, when the initial data are restricted to a subspace of the Hilbert 
space L 2 (il). We restrict ourselves to the weighted space 

L 2 (fl,K) with K(x) := e^/ 4 , (5) 

which means that the initial data are required to be sufficiently rapidly decaying 
at the infinity of the strip. As a measure of the additional decay, we consider 
the (polynomial) decay rate 

le :=sup{ 7 | 3C 7 >0,Vt>0, \\e {Ae+El)t \\ L2(nK) ^ L2{n) <C 7 (l+i)" 7 } • 

(6) 

Our main result reads as follows: 

Theorem 1. We have 70 = 1/4, while 7^ > 3/4. 

Here the power 1/4 corresponding to the untwisted case 9 = reflects the 
quasi-one-dimensional nature of our model (recall that c?/4 is the analogous 
decay rate for the heat semigroup in the L 2 -space over the whole Euclidean 
space R d for initial data in i 1 (M d ) D L 2 (R d , e^^dx)). Indeed, the result 
for 9 = follows easily by separation of variables (cf Section [3|). 

The essential content of Theorem [T] is that solutions to (PJ when the strip is 
twisted (i.e. 9 = tt) gain a further decay rate 1/2. The proof of this statement is 
more involved and constitutes the main body of the paper (cf Section [4]). It is 
based on the method of self-similar solutions developed in the whole Euclidean 
space by Escobedo and Kavian [5] and adapted to waveguide systems by the 
present authors in [TU], where it was shown that the heat kernel decays faster 
in geometrically twisted tubes than in untwisted ones. An open problem is to 
show that the decay rate 7^ is precisely 3/4 (cf Section [5]). 

The way how to understand the difference in the decay rates of Theorem [1] 
is due to a fine difference between the operators — Ao and —A n in the spectral 
setting: Although the operators have the same spectrum (cf ([3])), the shifted 
operator — Ao — E\ is critical, while — A ff — E\ is subcritical. The latter is 
reflected in the existence of a Hardy-type inequality 

- A OT - E x > p (7) 

with a positive function p (while such an inequality cannot hold for — Ao — E{). 

Various Hardy inequalities for — A„- — E\ were established in [8]. A general 
conjecture on the influence of the subcriticality of an operator on the improve- 
ment of the decay of the associated semigroup was made in jlOj , where an analog 
of Theorem[T]was proved for the decay rate in three-dimensional Dirichlet tubes. 
We also refer to [5] where the conjecture (for not necessarily self-adjoint opera- 
tors) is analysed from the point of view of heat kernels and its relationship with 
Davies' conjecture [2] is observed. 

The organization of this paper is as follows. In the following Section [2] we 
give a precise definition of the Laplacians — Ag and the associated semigroups. 
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The untwisted case is briefly treated in Section [3J obtaining, inter alia, the first 
statement of Theorem[T] The main body of the paper is represented by Section|4] 
where we develop the method of self-similar solutions to get the improved decay 
rate of Theorem [T] (and also to establish an alternative result, Theorem^. The 
paper is concluded in Section [5] by referring to physical interpretations of the 
result and to some open problems. 

2 Preliminaries 

The Laplacians — Ag are introduced as the self- adjoint operators associated on 
L 2 (fl) with the quadratic form ip H> ||V-0|| 2 having the domains 

2> 9 (n) := {ip e H\n) \ ip = o on rf } . 

Here and in the sequel || • || denotes the norm in the Hilbert space L 2 (Vt). It 
is possible to specify the operator domains (c/ [3]), but we will not need them. 
We only mention the result from [3J that the set of restrictions of functions from 
C^°(R 2 ) to ft that vanish on Tf is dense in Dg with respect to the 1 (£1) norm 
(cf m Ap P . B]). 

In view of ([3]), both the operators — Ag satisfy the Poincare-type inequality 

- A e > E x (8) 

in the sense of forms on L 2 (fl). Here E\ is the first eigenvalue of the one- 
dimensional operator — A^ D ^' a \ i.e. the Laplacian in L 2 ((— a, a)) subject to the 
Dirichlet boundary condition at —a and Neumann boundary condition at a. 
This inequality is sharp for — Ao, while it follows from [8] that (j8]) can be 
improved to a Hardy- type inequality ([7]) if 6 = tt. 

Recalling ((4]) and that we are interested in additional decay properties of ([T]), 
it is natural to rather consider the shifted parabolic problem (obtained from the 
standard heat equation (fTJ by the replacement u(x, t) M> e~ Elt u(x, t)): 

Ut-Au- Eiu = , (9) 

subject to the initial condition ([2]) and the boundary conditions 

u = on rfx(0,oo), ^ = ° on Te'xC ; 00 ); ( 10 ) 

where n denotes the normal vector to the boundary dfl. 

As usual, we consider the weak formulation of the problem and, with an 
abuse of notation, we denote by the same symbol u both the function on 
fl x (0, oo) and the mapping (0, oo) — > L 2 (£l). Standard semigroup theory 
implies that there exists a unique solution of (|9|- ([T0l) . subject to the initial 
condition ([2]), that belongs to C°([0,oo),£ 2 (f2)). More precisely, the solution 
is given by u(t) = Sg(t)uo, where 

Sg(t) := e ( A «+^)* (11) 
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is the heat semigroup associated with the shifted Laplacian — Ag — E\. 

It is easy to see that the real and imaginary parts of the solution u of ([9]) 
evolve separately. By writing u = + and solving © with initial 

data 9?(iio) and ^s(uq), we may therefore reduce the problem to the case of a 
real function uq, without restriction. Consequently, all the functional spaces are 
considered to be real in the sequel. 

3 The untwisted strip 

If the strip is untwisted (i.e. 9 = 0), the heat equation © can be easily solved 
by separation of variables. Indeed, the Laplacian — Ao can be identified with 
the decomposed operator 

(-A'j^l + l^f-A^) in L 2 (R)®L 2 ((-a,a)), (12) 

where — A R denotes the one-dimensional free Hamiltonian (i.e. the usual self- 
adjoint realization of the Laplacian in L 2 (R)) and 1 stands for the identity 
operators in the appropriate spaces. 

The eigenvalues and (normalized) eigenfunctions of — A^,^' a ' ) are respectively 
given by (n = 1,2,...) 

E n := (2n - l) 2 £i , J n (y 2 ) := ^ sin [E n (y 2 + a)} , (13) 

while the spectral resolution of — A R is obtained by the Fourier transform. Then 
it is easy to see that the heat semigroup So(t) is an integral operator with kernel 

oo 

s (x,x / ,t) :=Y,e~ iE ^ El)t Jn(x 2 )p(x 1 ,x' 1 ,t)J n (x / 2 ), (14) 

71=1 

where 

e -(a;i- a; ' 1 ) 2 /(4t) 

p(x 1 ,x 1 ,t) := ■= 

V47rf 

is the well known heat kernel of — A R . 

Using the explicit form of the heat kernel, it is straightforward to establish 
the following bounds: 

Proposition 1. There exists a constant C such that for every t > 1, 
£r-i f -i/4 ^ ||5 (t)|| L2(a!K) ^ i2 ( O ) < Ci -1 / 4 . 

Proof. To get the lower bound, we may restrict to the class of initial data 
of the form u (x) = ip(x 1 )J 1 (x 2 ) with <p G L 2 (R, Ki), K^xi) := e x ^ 4 . Then it 
is easy to see from (fT4)) that 

\\So(t)\\L 2 (n,K)^L 2 (n) > ||-P(*)|U 2 (B,Xi)->-L 2 (M) ! 



5 



where P(t) is the heat semigroup of — A K for which the lower bound with i -1 / 4 
is well known (or can be easily established by taking ip = K^ a with any a > 1/2 
and evaluating the integrals with the kernel p explicitly). On the other hand, 
using the Schwarz inequality, we have 

\\S (t)u \\ 2 < \\u f K f s^x'tfKix'^dxdx' 



oo „ 

\ V e-HEn-E^t / p ( Xl , x > lit ) 2 ifi(xi)- 1 dx x dx\ 



n=l 

for every uq E L 2 (il,K). Here the sum can be estimated by a constant in- 
dependent of t > 1 and the integral (computable explicitly) is proportional 

tor 1 / 2 . □ 

Remark 1. It is clear from the proof that the bounds hold in less restrictive 
weighted spaces. Indeed, it is enough to have a corresponding result for the 
one-dimensional heat semigroup P(t). 

As a consequence of Proposition [U we get: 
Corollary 1. We have 70 = 1/4. 

Proof. The lower bound of Proposition [1] implies 70 < 1/4. The opposite in- 
equality follows from the upper bound and Q. □ 



4 The self-similarity transformation 

Our method to study the asymptotic behaviour of the heat equation (|9|) is to 
adapt the technique of self-similar solutions used in the case of the heat equation 
in the whole Euclidean space by Escobedo and Kavian [S] to the present prob- 
lem. Following [ID] , devoted to the analysis of the heat kernel in twisted tubes, 
we perform the self-similarity transformation in the first (longitudinal) space 
variable only, while keeping the other (transverse) space variable unchanged. 



4.1 An equivalent time-dependent problem 

More precisely, we consider a unitary transformation U on L 2 (il) which as- 
sociates to every solution u £ L 2 OC ((0, 00), dt; L 2 (Q,dx)) of (JUJ) a self-similar 
solution u = Uu in a new s-time weighted space L 2 OC ((0, 00), e s ds; L 2 (0, <iy)) 
via 

u(Vl, J/2, s) = e s/i u{e s/2 yi ,y 2 ,e s - 1) . 
The inverse change of variables is given by 

u(x!,X2,t) = (t + u({t + l)- 1/2 Xl , x 2 ,log(t + 1)) . 

When evolution is posed in that context, y = (2/1,2/2) plays the role of space 
variable and s is the new time. 
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It is easy to check that, in the new variables, the evolution is governed by 

y-B — \yi 9iu - dfu - e s d\u -E 1 e s u~jU = Q (15) 

subject to the same initial and boundary conditions as u in ([2]) and (fTU|) . re- 
spectively. 

Remark 2. Note that (fl5|) is a parabolic equation with s-time-dependent coef- 
ficients. The same occurs and has been previously analyzed in twisted three- 
dimensional tubes [TU] and for a convection-diffusion equation in the whole space 
but with a variable diffusion coefficient [3] . A careful analysis of the behaviour 
of the underlying elliptic operators as s tends to infinity leads to a sharp decay 
rate for its solutions. 

Since U acts as a unitary transformation on L 2 (f2), it preserves the space 
norm of solutions of © and (|15[) . i.e., 

Il«(*)ll = l|fi(*)l|. (16) 

This means that we can analyse the asymptotic time behaviour of the former 
by studying the latter. 

However, the natural space to study the evolution (fT5|) is not L 2 (£l) but 
rather the weighted space ((5]). Following the approach of [10] based on a theorem 
of J. L. Lions [TJ Thm. X.9] about weak solutions of parabolic equations with 
time-dependent coefficients, it can be shown that (|15p is well posed in the scale 
of Hilbert spaces 

3) e (fi, K) c L 2 (f>, K) c S) fl (n, K)* , 

with 

S e (0,if) := {UeH^njq | u = on if}, 
where H 1 ^, K) denotes the usual weighted Sobolev space. 



4.2 Reduction to a spectral problem 

Multiplying the equation (p~5]) by uK and integrating by parts (precisely this 
means that we use uK as a test function in a weak formulation of (|15j) ). we 
arrive at the identity 

~ s \\u{s)f K = -J s e [u{s)]. (17) 
Here || • \\k denotes the norm in ((5|) and 

Jl\u] := WdMl + e s \\d 2 uf K E X e s \\uf K \ \\uf K 

is a closed quadratic form with domain D(Jg) := Qe(Q, K) (independent of s). 
It remains to analyse the coercivity of J| . 



7 



More precisely, as usual for energy estimates, we replace the right hand side 
of (fT7|) by the spectral bound, valid for each fixed s <E [0, oo), 

VueD(J$), J$[u]>m(s)\\u\\ 2 k , (18) 

where fJ-$(s) denotes the lowest point in the spectrum of the self-adjoint oper- 
ator Tg associated in L 2 (£l,K) with J|. Then p7|) together with ([T8")) implies 
the exponential bound 

V*€[0,oo), ||fi(«)||jr< ||uo||jce-J'o*wW*. (19) 

In this way, the problem is reduced to a spectral analysis of the family of oper- 
ators {T 9 S } S > . 

4.3 Study of the spectral problem 

In order to investigate the operator Tg in L 2 (Q, K), we first map it into a 
unitarily equivalent operator Tg :— UTgU^ 1 in L 2 (il) via the unitary transform 

K 1/2 u. 

By definition, Tg is the self-adjoint operator associated in L 2 (Vl) with the 
quadratic form Jg[v] := Jg{U~ 1 v], v £ D(Jg) := UD(Jg). A straightforward 
calculation yields 

J|M = \\d lV \\ 2 + i ||^|| 2 + e s ||d 2 t-|| 2 - B X e s |M| 2 , 

« e £>( J|) = s fl (n) n L 2 (r>, y \ dy) . 

In particular, D(Jg) is independent of s. In the distributional sense, we can 
write 

f$ = -d\ + y 2 - e s d 2 - B X e s . (21) 

We observe that the 'longitudinal part' of Tg coincides with the quantum 
harmonic-oscillator Hamiltonian 

H:=-^ + ^vl in L 2 (R) (22) 

(i.e. the Friedrichs extension of this operator initially defined on Cq°(M)). We 
recall the well known fact that the form domain 

D{H 1 ' 2 ) = H 1 (R)nL 2 (R,y 2 1 dy 1 ) 

is compactly embedded in L 2 (R), so that the spectrum of H is purely discrete. 
In fact, the spectrum can be computed explicitly (see any textbook on quantum 
mechanics, e.g., [7J Sec. 2.3]): 

<H) = \\(n+^))°° . (23) 

n=0 



2 V 2 
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Using now the discreteness of spectra of H and —Ajj^ ,a ' together with the 
minimax principle, one may easily conclude that also T| (and therefore Tjj) is 
an operator with compact resolvent for all s € [0, oo). In particular, ij.q(s) rep- 
resents the lowest eigenvalue of Tg . 

4.4 The asymptotic behaviour of the spectrum 

In order to study the decay rate via (fl~9|) . we need information about the limit 
of the eigenvalue /Ue(s) as the time s tends to infinity. Notice that the scaling 
of the transverse variable in (f21j) corresponds to considering the operator Tg 1 in 
the shrinking strip R x (— e~ s / 2 a, e~ s / 2 a). This suggests that Tg will converge, 
in a suitable sense, to a one-dimensional operator of the type (|22|) . We shall 
see that the difference between the twisted (0 = n) and untwisted case (9 = 0) 
consists in that the limit operator for the former is subject to an extra Dirichlet 
boundary condition at yi = 0. 

Thus, simultaneously to H introduced in (j2"2")l , let us therefore consider the 
self-adjoint operator Hp in L 2 (M) whose quadratic form acts in the same way 
as that of H but has a smaller domain 

D{H]> 2 ) := { V e DiH 1 ' 2 ) | p(0) = 0} . 

In fact, it is readily seen that Tq can be identified with the decomposed 
operator 

H®l + l®{-e s A { ~£ a) -E x e s ) in L 2 {R) <g> L 2 {{-a, a)) , (24) 

where 1 denotes the identity operators in the appropriate spaces. Using f|23[) . it 
follows that /Xo(s) = 1/4 for all s e [0, oo). Consequently, 

Aio(oo) := lim no(s) = 1/4. (25) 

s— >-oo 

Moreover, (|24l) can be used to show that TJ converges to H in the norm-resolvent 
sense as s — » oo, if the latter is considered as an operator acting on the subspace 
of L 2 (fi) consisting of functions of the form <p(yi)Ji(y%), where J\ is introduced 
in (jT3]l, 

It is more difficult (and more interesting) to establish the asymptotic be- 
haviour of fi n (s). A fine analysis of its behaviour leads to the key observation 
of the paper, ensuring a gain of 1/2 in the decay rate in the twisted case. 

We decompose the Hilbert space L 2 (f2) into an orthogonal sum 

L 2 (n) ^f)!®^, (26) 
where the subspace Sji consists of functions of the form 

^i(y) = v(»i)^i(8gn(-»i)»a). (27) 

Notice that 7/2 H> Jxi~ J/2) is an eigenfunction of — A^ N ^ a \ i.e. the Lapla- 
cian in L 2 ((— a, a)) subject to the Neumann boundary condition at —a and 
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Dirichlct boundary condition at a (reversed boundary conditions with respect 
to — A^ J ^' a ' ) ). Hence ipi satisfies the boundary conditions of — A T . Given any 
ip £ L 2 (51), we have the decomposition ip = ipx + with ipi £ fji as above and 
4> £ . The mapping i : ip H ► ^i is an isomorphism of i 2 (R) onto Hence, 
with an abuse of notations, we may identify any operator h on L 2 (K) with the 
operator uhi' 1 acting on fji C L 2 (il). Having this convention in mind, we state 
the following convergence result. 

Proposition 2. The operator T£ converges to Hp (BO 1 " in the strong-resolvent 
sense as s-> oo, i.e., 



VFe L 2 (n), 



lim 



i 



F ■ 



H 



D 



1 



©0 J 



F 



0. 



Proof. We proceed as in the proof of [TUJ Prop. 5.4]. For any fixed F £ L 2 (J7), 
let us set ip s := (T£ + In other words, ip s satisfies the resolvent equation 



Vv£D(j°), j*(v,iP s ) + (v,iP s ) = (v,F), 



(28) 



where (•,•) denotes the inner product in L 2 (il) and J*(-,-) is the sesquilinear 
form associated with (|20j) . In particular, choosing ip s for the test function v 
in (|2"51). we have 



WW. 



+ jglltfi^ll 5 



£i||^ 



W^^ll^f + M 2 - (29) 



Notice that 1 1 1 1 2 > -E'lllV'sll 2 by Fubini's theorem and the Poincare inequality 



for - A 



( — a. a) 
DN 



and — A 



( — a, a) 



. Consequently, 



ND 

\\diip s f<C, \\yiA\\ 2 <C, Us\\ 2 <C, ||ft^|| 2 -Bi||^|| 

where C is a constant proportional to ||-F|| 2 . 
Now we employ the decomposition 



V'a(y) = <fs(yi) Ji(sgn(-?/i)?/2) +</> s (y) 



where </> s £ fj^, i.e., 

Vj/i e k 



Ji(sgn(-yi)2/ 2 ) <t>s(yi,y2)dy 2 = 0. 



(30) 



(31) 



That is, y2 s (2/1,2/2) is orthogonal to the ground-state eigenfunction of 



-A 



( — a, a) 



D jy (respectively of - A^ 



(-0,0)n 



if j/i < (respectively j/i > 0). Then 



|d 2 V s || 2 -£i||V s || 2 = ||». 



2<Ps\ 



12 1| |^ji 2 -.e 

> \ \\d 2 (jys\\ 2 + {\E 2 -E 1 



2 1 { 1 TP \ II A. II 2 
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where E2 = 9Ei denotes the second eigenvalue of — A^,^'^ (which coincides 
with that of —A N p ). Thus it follows from the last inequality of (|30|) that 

\\<t> s \\ 2 < Ce~ s and \\d 2 (f>s\\ 2 < Ce' 3 , (32) 

where C is a constant proportional to ||F || 2 . 

It follows from (|30|) that {ip s }s>o is a bounded family in D(J 3 ). Therefore it 
is precompact in the weak topology of D(J 3 ). Let "0oo be a weak limit point, i.e., 
for an increasing sequence of positive numbers {s„}„eN such that s n — > 00 as 
n — > 00, {i/j Sn } n £Fi converges weakly to ipoa in D(J 3 ). Actually, we may assume 
that it converges strongly in L 2 (Vl) because D(J 3 ) is compactly embedded in 
L 2 (fl). Since {(j) Sn } n ^ converges strongly to zero in L 2 (f2) due to (|32| . we 
know that t/j^ £ i.e., 

ipoo(y) = V>oo(yi)i7i(sgn(-2/i)y 2 ) 

with some ifioo G L 2 (M.). Since the weak derivative di'f/'oo € L 2 (il) exists, we 
necessarily have ip oo € if 1 (R) and 

<^oo(0) =0. 

Finally, let ip € C£°(R\{0}) be arbitrary. Taking 

v(y) := ip(y 1 )J 1 (sgn(-y 1 )y 2 ) 

as the test function in (|28| . with s being replaced by s„, and sending n to 
infinity, we easily check that 

(<P, </?oo)l 2 (R) + — (yxLp, yiipoo)h^(WL) + {<P, ^oo)L2( R ) = (ip, f) L ^(M) , 

where a 

f(yi)-= / Ji(sgn(-yi)y 2 ) F{ Vl ,y 2 )dy 2 . 

J —a 

That is, (foo = (Hp + l) -1 /, for any weak limit point of {tp s } s >Q. Summing 
up, we have shown that ip s converges strongly to ipoa in L 2 (Q) as s — > 00, where 
il> oo (v)=[(H D +z)- 1 ®0 ± ]F. □ 

Corollary 2. One has 

/i w (oo) := lim /i(s) =3/4. 

5—)- OO 

Proof. In general, the strong-resolvent convergence of Proposition [2] is not suf- 
ficient to guarantee the convergence of spectra. However, in our case, since the 
spectra are purely discrete, the eigenprojections converge even in norm (cf \12\). 
In particular, (s) converges to the first eigenvalue of Hd as s — > 00. It remains 
to notice that the first eigenvalue of Hd coincides (in view of the symmetry) 
with the second eigenvalue of H which is 3/4 due to (j2"3")l . □ 



11 



4.5 A lower bound to the decay rate 

We come back to (fTO]) . Recalling ([21)1 and Corollary [3J we know that for arbi- 
trarily small positive number e there exists a (large) positive time s £ such that 
for all s > s £ , we have /ig(s) > /ie(oo) — e. Hence, fixing e > 0, for all s > s e , 
we have 



fi 6 (r) dr < - (in(r) dr-[pig(oo) - e](s - s e ) 
o Jo 

< \pe(po) - e]s e -[^ e (oo) - e]s, 

where the second inequality is due to the fact that ^e(s) is non- negative for all 
s > (it is in fact greater than or equal to 1/4, cf Proposition [3] below). At the 
same time, assuming e < 1/4, we trivially have 

- / fig(r) dr < < [/ie(oo) - e]s £ ~[fig(oo) - e]s 



o 



also for all s < s e . Summing up, for every s £ [0,oo), we have 

\\u(s)\\ K <C e e-^°°^ s \\u \\ K , (33) 

where C e := e s ' > e \»i>(°°)-e]<>e , 

Now we return to the original variables (xi,X2,t) — (e s / 2 yi,y2,e s — 1). 
Using (|T51) together with the point-wise estimate 1 < K, and recalling that 
uq = i*o, it follows from (|33[) that 



||«(*)ll = < l|fi(*)lk <Ce(l + i)-^(°°)- £ ] || Uo || A - 

for every t £ [0, oo). Consequently, we conclude with 

\\S e (t)\\L^n,K)^mn) = sup < C £ (1 + 1 )-^(<»)-=] 

UDei 2 (n,if)\{o} \\ u o\\k 

for every t £ [0,oo). Since e can be made arbitrarily small, this bound implies 

le > Me(oo) . 

This together with Corollary [T] proves Theorem [TJ 

4.6 A global upper bound to the heat semigroup 

Theorem [T] provides quite precise information about the extra polynomial decay 
of solutions u of © in a twisted tube in the sense that the decay rate 7^ is 
better by a factor 1 /2 than in the untwisted case. On the other hand, we have no 
control over the constant C 7 in © (in principle it may blow up as 7 — > j g ) . We 
therefore conclude this section by establishing a global (in time) upper bound 
to the heat semigroup (i.e. we get rid of the constant C 7 ) but the prize we pay 
is just a qualitative knowledge about the decay rate. It is a consequence of (fT9t 
and the following result: 
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Proposition 3. Vs > 0, Mo(s) = 1/4, fi„(s) > 1/4. 

Proof. The identity for fi is readily seen from the decomposition and (|2U)) . 
Using Fubini's theorem and the minimax principle, it is also easy to deduce 
from (|2T))) that fJ, n (s) > 1/4 for all s > 0. To show that the inequality is strict, 
let us assume by contradiction that /^(s) = 1/4 for some s > 0. Let v denote 
the corresponding eigenfunction of T*. Then the identity = /i7r(s)||u|| 2 

yields 

IIMI a + ^IM| 2 = 7W a and \\d 2 vf-E 1 \\v\\ 2 =0. (34) 
16 4 

Using the direct-sum decomposition (|26|l . the second identity implies that w is 
of the form (|27[) . The continuity of the eigenfunction i> inside f2 in turn requires 
that ip(0) = 0. However, this contradicts the first identity in (|3"4"|) which says, 
in view of (|23[) . that is the first (therefore nowhere vanishing) eigenfunction 
of H. □ 

Combining this result with Corollary [21 we see that the number 

c e := inf /ifl(s)-l/4 (35) 

sS[0,oo) 

is positive if 9 = tt and zero if 9 = 0. In any case, (|19l) implies 

||fi(«)||jr < \\uo\\i<e-^ +1 ^ s 

for every s£ [0, oo). Using this estimate instead of ff33|) . but following the same 
type of arguments as in Section 14.51 below ([33| , we thus conclude with: 

Theorem 2. We have 

Vt e [0, co), \\S e (t)\\ L 2 (njK) ^ L2{n) < (l + i)- (C8+1/4) , 

where > ("and Co = 0J. 

5 Conclusions 

The classical interpretation of the heat equation (p} is that its solution u gives 
the evolution of the temperature distribution of a medium in the strip 51 sur- 
rounded by a perfect insulator on the Neumann boundary and by a substance 
of a high- heat capacity and of zero temperature on the Dirichlet boundary . It 
also represents the simplest version of the stochastic Fokker-Planck equation de- 
scribing the Brownian motion in ft which is normally reflected on and killed 
on (c/ [TT| for a probabilistic setting in an analogous higher-dimensional 
model). Then the results of the present paper can be interpreted as that the 
twisting of boundary conditions (i.e. 9 = ir) implies a faster cool-down/death 
of the medium/Brownian particle in the strip. Many other diffusive processes 
in nature are governed by (p}. 
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Our proof that there is an extra decay rate for solutions of ([T]) if the boundary 
conditions arc twisted was far from being straightforward. This is a bit surpris- 
ing because the result is quite expectable from the physical interpretation, if one 
notices that the twist makes it more difficult for the Brownian particle to pass 
through the channel at {x% = 0}, because of the proximity of killing boundary 
conditions. At the same time, the Hardy inequality ([7]) did not play any role in 
the proof of Theorem[T] (although, combining the theorem with the results of [8] , 
we eventually know that the existence of the Hardy inequality is equivalent to 
the extra decay rate for the heat semigroup). It would be desirable to find a 
more direct proof of Theorem Q] based on ([7]). 

We conjecture that the inequality of Theorem Q] can be replaced by equality, 
i.e., 7^ = 3/4 for twisted strip. The question of optimal value of the constant c v 
(and its quantitative dependence on the half-width a) from Theorem [2] also 
constitutes an interesting open problem. Note that the two quantities are related 
by + 1/4 < 7e . 

The present paper can be viewed as a continuation of the research initiated 
by our work [10] . where we investigated the large-time behaviour of the heat 
semigroup in geometrically twisted Dirichlet tubes. It confirms that the effect 
of twisting (leading to the subcriticality of the Laplacian and implying an im- 
provement of the decay rate of the associated heat semigroup) is more general, 
namely it holds true also in waveguide systems twisted via boundary conditions. 
We expect that the extra decay rate will be induced also in the systems twisted 
via embedding of the strip into a negatively curved manifold, for which the 
existence of Hardy inequalities is already known [9] . 

More generally, recall that we expect that there is always an improvement 
of the decay rate for the heat semigroup of a subcritical operator (c/ [TUJ Con- 
jecture in Sec. 6] and [51 Conjecture 1]). 
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